Black/ White hole radiation from dispersive theories 
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We study the fluxes emitted by black holes when using dispersive field theories. We work with 
stationary one dimensional backgrounds which are asymptotically fiat on both sides of the horizon. 
The asymptotic fiuxes are governed by a 3 x 3 Bogoliubov transformation. The fluxes emitted by the 
corresponding white holes are regular and governed by the inverse transformation. We numerically 
compute the spectral properties of these fluxes for both sub- and superluminal quartic dispersion. 
The leading deviations with respect to the dispersionless flux are computed and shown to be governed 
by a critical frequency above which there is no radiation. Unlike the UV scale governing dispersion, 
its value critically depends on the asymptotic properties of the background. We also study the 
flux outside the robust regime. In particular we show that its low-frequency part remains almost 
thermal but with a temperature which significantly differs from the standard one. Application to 
four dimensional black holes and Bose Einstein condensates are in preparation. 
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I. INTRODUCTION 

In 1981, Unruh showed [1] that the phonon field equa- 
tion in a nonhomogeneous flow is analogous to that of a 
relativistic field propagating in a four dimensional curved 
space-time. He also suggested that this analogy could 
be used to experiment black hole (BH) radiation in the 
lab. However, as pointed out in [2], this analogy is lim- 
ited by the fact that in condensed matter the phonon 
field becomes dispersive above a certain wave vector, say 
A. Motivated by this remark, Unruh numerically showed 
that the fluxes emitted by an acoustic black hole are not 
significantly affected by the dispersive properties of the 
phonons, at least when A 3> k, where k, is the gradient 
of the fluid velocity at the sonic horizon. [3] Since then, 
this robustness of the fluxes has been confirmed both by 
algebraic [4-9] and numerical [10, II] techniques. 

However, many open questions remain. In particular, 
little is known about 

• the exact nature and the precise range of the pa- 
rameters delimiting the robust regime, 

• the scaling properties of the leading deviations in 
the robust regime, 

• the properties of the fluxes outside this regime, and 

• the comparison of sub- and superluminal disper- 
sion. 

In view of the possibility to detect Hawking radiation in 
the lab, [12] it is important to provide quantitative esti- 
mates of the modifications of the fluxes induced by dis- 
persion. Besides this reason, answering the above ques- 
tions is also relevant when exploring the consequences 
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on black hole physics due to violations of Lorentz invari- 
ance [2, 13-16] that could stem from (or be related to) 
quantum gravity. This could be true in particular when 
considering the signals emitted during the evaporation of 
light black holes (M ~ 10^^^ Planck mass) that could be 
produced at the LHC [17]. 



In this paper, we aim to answer these questions by a 
two step analysis. First we provide a complete descrip- 
tion of the quantum propagation of a dispersive field in 
stationary black hole geometries. (To our knowledge, this 
has not been presented anywhere.) Secondly, focussing 
on quartic dispersion and using numerical methods, we 
study both the small modifications of the fluxes induced 
by dispersion when A 3> k, and the larger modifications 
when this inequality is not satisfied. 



We are planning to extend this work in two directions. 
First, to Bose condensates, where the phonon field obeys 
the Bogoliubov-de-Gennes equation which differs from 
the "scalar" field equation we consider below (this pro- 
gram has been completed during the publication process 
of the present work [18]) and secondly to 4D black holes, 
where the gravitational potential induces grey body fac- 
tors which could play a critical role. 



In Sec. II, we present the class of velocity profiles and 
the stationary mode equation. In Sec. Ill, we study the 
algebraic properties of the modes and the Bogoliubov 
transformation one is dealing with. In Sec. IV, we present 
our numerical treatment. In Sees. V and VI, we numer- 
ically solve the mode equation and interpret the results 
for sub- and superluminal quartic dispersion respectively. 
We conclude in Sec. VII. 



II. THE SETTINGS 



A. Velocity profiles 



We work in two space-time dimensions and consider 
stationary black hole geometries. These shall be charac- 
terized by the velocity profile v{x): 



ds'^ = -di^ + {dx - v{x)dtf 



(1) 



This expression appears [1] when considering the propa- 
gation of low-frequency phonons in a moving fluid whose 
velocity field v{x) is measured in the Galilean frame (not 
necessarily the lab frame) where v only depends on x. 
The sound velocity is assumed to be constant and has 
been set to 1. A more general case will be considered 
in [18]. The above metric possesses an event horizon 
where the flow becomes supersonic. 

When the fluid flows to the left, ?; < 0, and when 
\v\ increases toward the left, one obtains a future (black 
hole BH) horizon. Had we considered the flow —v{x), 
we would have obtained a white hole (WH) geometry. 
As we shall see, the spectral properties (in terms of the 
conserved frequency lo = idt) of fluxes emitted by this 
WH are in one-to-one correspondence with those of the 
fluxes emitted by the BH governed by v{x). However 
they differ slightly. It should also be pointed out that, 
since we only consider dispersive mode; equations, sec e.g. 
Eq. (5), — oo < t,a; < oo in Eq. (1) represents the whole 
space-time [19], and not only a part of it as would have 
been the case when using relativistic fields. 

We shall also assume that v becomes asymptotically 
constant on both sides of the horizon. The existence of 
two asymptotic regions is important because one then has 
well-defined asymptotic modes on both sides, and there- 
fore well-defined particle fluxes. This is particularly rel- 
evant for the proposal [20, 21] based on measurements of 
the long distance correlations [22, 23] between Hawking 
quanta and their partners propagating on the other side 
of the horizon. 

In order to determine the generic consequences of dis- 
persion, we shall consider a large class of flows. In con- 
trast with the nondispersive case where they play no role, 
we shall see that the asymptotic properties of v{x) are 
crucial when using nonlinear dispersion relations. This 
will imply in particular that the inequality A 3> At is not 
sufficient to guarantee no significant deviation from the 
standard fluxes. 

The class we shall use generalizes that of [10] and con- 
tains profiles of the form 



v{x) = -1 + Ds\gn{x) tanh^/" 



D 



(2) 



When changing n and D, k, the slope of v at the hori- 
zon, is fixed, as is the location of the horizon at x = 0. 
When using the 2D massless relativistic theory, the emit- 
ted spectrum only depends on n. Instead, when using 
dispersive theories, together with the ratio A/k, both n 




Figure 1: Velocity profile 11(3;; /t, n, D) as a function of kx. (a) 
n fixed to 2, D = 0.5 (solid line) and D = 1 (dashed line), 
(b) D fixed to 0.5, n = 1 (solid line), n = 20 (dashed line). 



and D also affect the properties of the fluxes. Figure 1 
shows v{x) for several values of D and n, as a function 

of KX. 

The power n controls the sharpness of the transition 
from the linear behavior of v near the horizon to the 
asymptotic flat region. For n ^ 00, the transition be- 
comes sharp. When using dispersive theories, sharp tran- 
sitions give rise to nonadiabatic effects which produce su- 
perimposed oscillations [10]. Unless said otherwise, n is 
equal to 2 in the following. 

The parameter D e]0, 1] flxes the asymptotic values 
of the velocities, 



v± 



-1±D, 



(3) 



for a; —> ±00 respectively. When using nonlinear disper- 
sion relations, D flxes the critical frequency Wmax above 
which the radiation vanishes. This is relevant because D 
is constrained to be much smaller than 1 in some propos- 
als of experiments [24] . More generally, we shall flnd that 
t^max is the most relevant parameter, in that the correc- 
tions induced by dispersion are essentially governed by 
K/wmax, and not by k/A as one might have expected. 



B. Dispersive wave equation 

In Nature of course, various dispersion relations are 
found depending on the condensed material used. How- 
ever, when imposing analyticity in fc^, quartic dispersion 
relations represent the first two possible deviations from 
the linear case. (Dissipative theories constitute another 
class that should be considered separately [25, 26].) In 
this paper, we shall thus restrict our numerical analysis 
to 



02 = fc2 ± 



A2' 



(4) 
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where ft and k are the energy and momentum measured 
in the frame comoving with the flow. 

When working with the metric of Eq. (1) the cor- 
responding wave equation for the velocity potential (j) 
reads [3] 

{dt + d,v){dt + = {dl T ^dt)<t>. (5) 

when assuming that the fluid density is constant. (This 
restriction will also be removed in [18].) The upper 
(lower) sign corresponds to super- (sub-) luminal disper- 
sion. The conserved scalar product takes the form 

/oo 
dx [4>l{dt + vdx)<i)2 - hidt + vda:)(l)l] . 
-oo 

(6) 

It is independent of the dispersion and it coincides with 
the standard Klein-Gordon (KG) product evaluated on 
the preferred time slices specified by dispersion (here t = 
est). 

Since the flow is stationary, we shall work at fixed 

(Killing) frequency oj. Setting (j) = e~^*<j)u{x) in the 
wave equation yields 

{-iw + dxv){-iw + vdx)(pcj = [dl T ■^i9;^)'/>a;- (7) 

As usual, the scalar product is block diagonal in oj. 
Therefore, when working with stationary states, all ob- 
servables will be expressed as sums of observables defined 
at fixed ijj. More precisely, one should consider the couple 
of frequencies (w, — w) since opposite values of w mix in 
stationary Bogoliubov transformations. In what follows, 
to avoid any confusion, lo is always positive definite. 

Before numerically studying the solutions of Eq. (7), it 
is worth analyzing the enlarged set of solutions and the 
Bogoliubov transformation one obtains. 

We remind the reader that the oscillatory solutions to 
Eq. (7) can be classified according to the sign of their 
group velocity in the comoving frame, = dCl/dk. In 
the following we call the modes with Vg""^ < left movers 
and those with the opposite sign right movers. In the ab- 
sence of dispersion (and for 2D massless fields), these two 
sectors are decoupled, and the Bogoliubov transforma- 
tion is particularly simple. When introducing (quartic) 
dispersion, four independent solutions to Eq. (7) exist. 
The classification between left and right movers can still 
be done, but the left-right decoupling is lost. It is worth 
mentioning that, when considering a varying speed of 
sound rather than a varying flow- velocity, see e.g. [21], 
this decoupling is lost even in the long wavelength ap- 
proximation {i.e. without dispersion). Notice also that 
one can introduce dispersion in a way that preserves this 
factorization [4, 27], but these models do not seem to 
govern condensed matter systems. 
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Figure 2: The two straight lines represent lo — V-yk for the 
subsonic one and uj — v-k for the supersonic one. The roots 
of Eq. (8) correspond to the abscissa of the intersections of 
these lines with the curves ±n(fc). The numerical values are 
p = 1 (quartic dispersion), A/k = 50, D = 0.9, w/zt = 10. 

III. ANALYTIC STUDY 

In several aspects, the material presented below closely 
follows what has been presented in [10]. The novelties are 
related to the fact that we treated both asymptotic re- 
gions on an equal footing. In particular, we establish the 
complete character of the "in" and "out" mode bases, and 
we write the unitary Bogoliubov transformation which 
governs the general case (for stationary, one dimensional, 
single horizon geometries). Finally, we show how to re- 
late the fluxes emitted by black and white holes. 

A. Asymptotic solutions 

Given our velocity profile, v is asymptotically constant 
in both regions \kx\ D. There, the solutions of Eq. (7) 
are superpositions of plane waves e**^^ with constant am- 
plitudes. To characterize a solution, one needs the roots 
k(ijj) of the asymptotic wave vectors and the amplitudes 
Ak of the corresponding waves. Owing to the nonlin- 
earity of the dispersion relation, in addition to the os- 
cillatory modes, exponentially growing/decaying modes 
exist, governed by complex values of k. 

We first study how to handle these extra modes in the 
present settings, paying attention to the completeness of 
the mode basis. To be specific, we concentrate on super- 
luminal cases, and only indicate the (small) differences 
that arise in subluminal cases. Let us thus analyze the 
roots of 

ju2p-|-2 

{^-v^kf = e + ^ = VL\k). (8) 

for a; > 0, and for an arbitrary value of the integer p > 0. 
We work with a polynomial form so as to know the num- 
ber of complex roots of Eq. (8). For everything concern- 
ing the real roots, nothing changes if one uses any strictly 
convex function f2^(fc) with a slope equal to 1 for low fc. 
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In the subsonie region, < 1, as can be seen in 
Fig. 2, two real roots exist: fc" > and fc^ < which cor- 
respond to a right and a left mover respectively. (When- 
ever an ambiguity could exist, we shall add a superscript 
u or to designate right or left movers.) p pairs of com- 
plex conjugated roots also exist, since the equation is 
real. One thus has p roots with a positive imaginary 
part (which correspond to growing modes to the right) 
and p solutions with a negative imaginary part (decaying 
modes) . 

In the supersonic region > 1, for w smaller than 

a critical frequency Wmax (we shall compute its value be- 
low) there exist four real roots (see Fig. 2), and thus only 
p — 1 pairs of complex solutions. In other words, when 
comparing this with what prevailed in the subsonic re- 
gion, a pair of growing and decaying modes has been 
replaced by a couple of oscillatory modes. Such replace- 
ments of pairs of modes is a generic feature of QFT in 
external fields, see the Appendix of [28]. In the present 
hydrodynamical settings it will occur at all horizons (for 
both black and white holes). For flows to the left, v <Q, 
the two new real roots both correspond to right movers. 
Notice that for subluminal dispersion, the replacement 
occurs the other way around since the four real roots are 
found where the flow is subsonic. 

Now, given that there are 2p + 2 roots, the general 
solution (t)uj{x) is a superposition of 2p-|- 2 modes. How- 
ever, when considering the Fourier transform of the field 
operator, (j)oj{x), only a subclass of these should be used. 
And this subclass is complete in a precise sense we ex- 
plain below. (To our knowledge, this crucial aspect has 
not been discussed in the context of acoustic black holes.) 
To establish the completeness of the mode basis we first 
need to compute the frequency aJmax which is introduced 
by both sub- and superluminal dispersion, and which will 
cut off Hawking radiation. 

B. Maximal frequency 

The frequency Wmax is the value of uj where the two 
extra rc;al roots merge into each other. It is thus reached 
when the straight Ymeu) — V- k is tangent to —^{k), as il- 
lustrated in Fig. 3, or equally when — \u\—V-k\s tangent 
to J7, where Q, is the square root of the rhs of Eq. (8). 

Restricting attention to quartic dispersion, the corre- 
sponding value of k is 



k^,^,+ = ^^v-'_-A+\v.\^v-'_+S. (9) 

The frequency Wmax is then obtained from Eq. (8). It is 
thus of the form 

Wmax,+ = A/+(£>). (10) 

For £) < 1, one has /+(£>) oc D^/^. This means that, 

whatever the value of the dispersion scale A, tiimax,-!- can 
be arbitrarily small, and in particular it could be even 
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Figure 3: Solutions to Eq. (8) in the supersonic region, for 

different values of u. For uj < (j^max/K, two roots exist in the 
bottom left quadrant. They move closer to one another as u) 
increases, and merge when the straight lino u) — kv- becomes 
tangent to — f2(fe), for uj = Wmax- For ui > Wmax, these roots 
no longer exist. 

smaller than the temperature K;/27r one would get in a 
dispersionless theory. 

For subluminal dispersion, the maximal frequency 
t^max.- is again defined by the merging of the two ex- 
tra real roots. In this case it occurs where \v^\ < 1. In 
the quartic case the corresponding value of k is 



fcmax,- = ^^4-V%-\v+\^vl + 8. (11) 

Thus Wmax,- is also of the form Wmax,- = ^f-{D). For 

< 1, we again obtain f-{D) oc D^/'^. 

Fig. 4 shows the contours of constant dJmax in the plane 
(Z), A), for sub- and superluminal dispersion. The same 
value of Wmax/K can be obtained with very different com- 
binations of {D,A/k). As we shall see in Sec. V, the 
fluxes are highly degenerate on these contours, which 
means that the corrections with respect to the standard 
fluxes are essentially governed by Wmax/w- 

C. Mode orthonormality and mode completeness 

To proceed to the canonical quantization, one should 
use mode bases that are orthonormal and complete. 
When the background geometry is homogeneous, this is 
rather easy, even in the presence of dispersion. However, 
when V of Eq. (1) varies and in particular when it char- 
acterizes a horizon, it becomes trickier. Let us thus first 
consider situations where v is constant both in space and 
time. 



1. Homogeneous metrics, k -representation 

In this case one should exploit the homogeneity and 
express the field operator in terms of exponentials e**^^ 
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Figure 4: Contours of constant Wmax in the plane (-D, A/k). D 
is the horizontal coordinate, A/k the vertical one. Prom left to 
right, the contours correspond to Wmax/w = 0.5, 2, 5, 10. The 
solid contours are for superluminal dispersion and the dashed 
ones for subluminal dispersion. The contours coincide for 
small D, and so do the functions /+ and /_. Their difference 
becomes visible as D grows. 



and creation/destruction operators labelled by (real) k 



(t>{t,x) = 



f 

J — ( 



dk 



ake 



Ikx 



-iujt 



+ h.c. 



(12) 



where u) = Vt{k) + vk with Vi{k) given in Eq. (8). In this 
representation, creation/destruction operators obey the 
usual commutators [a^, a\,\ = 6{k — k'). In addition, the 
modes which multiply these operators in Eq. (12) are or- 
thonormal: their norm, Eq. (6), is S(k — k'), irrespectively 
of both whether v is sub- or supersonic, and the choice 
of the dispersion relation (8).^ 

The completeness of the mode basis is needed to verify 
that the equal time commutator 



[(f>{t, x),'jf{t, x')] = iS{x — x'), 



is satisfied when (j) is given in Eq. (12) and where tt = 
{dt — vdx)4>- In the present representation, the proof is 
easily obtained by making use of the completeness (in the 



^ In the case of subluminal dispersion, there eould exist a maximal 
value of k associated with the vanishing of f2^(fc). In this case, 
the mode basis is no longer complete in a Fourier sense. To deal 
with this, one can either work [29] with "regularized" dispersion 
relations which do not develop imaginary values of f2, or with a 
"Bloch waves" perspective [30] where O is periodic in fc-space. 



sense of Fourier analysis) of the exponentials e*'^^, with k 
real from [—00,00]. It should be stressed that the proof 
applies to all dispersion relations and to all (constant) 
values of v, both sub- and supersonic. 

When considering nonhomogeneous stationary flows 
one must use the conserved frequency lo in place of the 
wave vector k to label modes and operators. So, in prepa- 
ration for this case, let us study the change of represen- 
tation from fc to a; in homogeneous flows. 



2. Homogeneous metrics, the u> -representation 

The above conclusion concerning the completeness 

character of the modes implies that, in any homogeneous 
flow, the "extra" (growing and decaying) solutions of the 
mode equation at fixed oj must be discarded when quan- 
tizing (j). When the flow is subsonic, using the disper- 
sion relation (8), one should thus discard the 2p grow- 
ing/decaying modes, whereas when the flow is supersonic, 
one should discard only 2(p— 1) modes. Because of this, 
the change of representation from fc to a; should be done 
separately in sub- and supersonic regions. 

Let us first consider the subsonic case where only two 
real roots exist. For each w > 0, one has two modes: 




(14) 



where A;* = k'^{uj), i = u,v, are the two real roots 
of Eq. (8) describing the right- and left-moving modes. 
These modes obey Eq. (7) and are orthonormal for the 
KG product, in the sense of a Dirac distribution S{ui—ui'). 
Correspondingly, one has the rescaled annihilation oper- 
ators: 



(15) 



which verify [a^, d^t] = 6^^ d{uj — uj'). Using these modes, 
the integral in Eq. (12) can be rewritten as 




4){t,x) 



doj 



h.c. 



= ^l{x)dl + ^l{x)al 



(16) 
(17) 



(13) The family of modes (14) is complete, because the expo- 
nentials e^^^ are. 

When the flow is supersonic, the situation is trickier 
because the jacobian dk/dw crosses zero for right movers. 
Nevertheless, similar results are obtained. Namely, as 
above, one first separates the integral dk in Eq. (12) 
into two integrals over a u and a v sector. When consider- 
ing the V sector {k < 0) in left-moving flows v < 0, noth- 
ing changes because dk/dio does not cross zero. There- 
fore, all left-moving (positive norm) modes can still be 
uniquely described by cj belonging to [0, 00], and the 
change of variable k ^ ui can be made without further 
precaution. 
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The same is no longer true for the right-moving sector. 
In fact, the integral over the right movers dk splits 
into an integral over co belonging to [0, oo] plus another 
piece over negative frequencies belonging to [— Wmax,0]. 
In addition, for a given value of w < in this interval, 
two real roots > exist. It ends when the two 

roots merge into each other when — Wmax is reached. 

Thus, for u) > Wmax, one has only one (positive norm) 
u-root and ^u> reads as in Eq. (17). Instead, when < 
<^ < i^max, three real roots exist: the continuation (in iv) 
of the former positive norm one, and two new roots with 
negative O, and thus negative norm, see Fig. 3. In this 
case, 0a; must be decomposed as 

1=1,2 

When compared with Eq. (17), the last two terms de- 
scribe the two new roots. A complex conjugate and a 
subscript — have been used to characterize the modes 
which multiply the creation operators al^ This means 
that the modes ip'^^ ; have a positive norm and obey the 
mode equation with a frequency idt = -co < 0. 

Using the Fourier operator (j)u;{x) of Eq. (18), con- 
taining both annihilation and creation sectors, the field 
(j){t,x) can still be written as in Eq. (16), as an integral 
over w e [0, oo]. 

3. Inhomogeneous metrics 

In metrics which contain a transition from a subsonic 
to a supersonic flow, the decomposition of (j)^ is modified: 
because of the scattering on v{x), the asymptotic modes 
used above mix into each other, and with the growing 
and decaying modes. As explained in Appendix A, when 
excluding the growing modes, there is a change in the 
dimensionality of the mode bases. Indeed, instead of 
Eq. (18), for < w < Wmax, 0w now reads 

= al ^l{x) + al ^l{x) + ati {<pl^{x)Y . (19) 

We have introduced only one extra mode, and not two, 
because the combination of the former two (^"^ ; that is 
orthogonal to the above (p^^ in the asymptotic super- 
sonic region is unbounded in the subsonic region. This is 
very similar to what is found with the Airy function: only 
one combination of two oscillatory modes in the allowed 
region remains bounded in the forbidden region. 

For LO > Wmax iustcad, nothing changes and 4>i^ should 
still be decomposed as in Eq. (17) because one has two 
(positive norm) real roots for all x. 

When considering subluminal dispersion, the situation 
is similar but with several changes. For < lu < Wmax, 
there are four real roots in the asymptotic subsonic part 
of the flow, and the mode operator should be decom- 
posed again as in Eq. (19). The difference with super- 
luminal dispersion is that the two additional roots are 



associated with the first term a" and no longer with 
the third one. 



D. Bogoliubov transformation 

For stationary flows and stationary states, the proper- 
ties of the fluxes are all encoded in a (w-block-diagonal) 
Bogoliubov transformation. This relates the in modes 
j^rn^^ which are associated with the destruction op- 
erators annihilating the initial vacuum state, to the out 
modes , (f°"^ which characterize the asymptotic par- 
ticle content of the fluxes. 

To obtain the asymptotic temporal behavior of the 
modes from their asymptotic spatial behavior, it suf- 
fices (see e.g. section 1.3. in [23]) to consider wave- 
packets centered about the value of u under examina- 
tion. In what follows we shall not distinguish single 
frequency modes from the corresponding broad wave- 
packet. Therefore when discussing the initial (or final) 
behavior of a mode, it should be understood as that of 
the corresponding wave packet. 

Given that the flow v becomes asymptotically constant 
for large |a;|, the identiflcation of in and out modes is free 
of ambiguity. The in (resp. out) modes are deflned as the 
modes with positive comoving frequency CI in the asymp- 
totic past (resp. future) and in the regions where the flow 
is homogeneous. They are separated into u, in modes 
(resp. u,out modes) with dCt/dk > in the asymptotic 
past (resp. future), and v,in (resp. v,out) modes with 
dfl/dk < 0. These conditions fix uniquely the complete 
set of 3 in (resp. out) modes. In this work, the state of 
the field is taken to be the in- vacuum, annihilated by the 
a™ associated with the in-modes. Physically, it corre- 
sponds to a state devoid of particles propagating toward 
the horizon in the regions where v is constant. Notice 
finally that this definition of the in vacuum docs not co- 
incide with that of the "Unruh- vacuum" [31] which rests 
on the use of relativistic fields, and of an afline null pa- 
rameter on the past horizon. 

1. Simplified case 

To sec what the new aspects brought in by dispersive 
effects are, let us briefly describe the Bogoliubov transfor- 
mation when assuming that the left-moving piece a^0^ 
in Eqs. (17, 19) decouples from the right-moving sector. 
In this case it is sufficient to consider only one relation 
amongst u modes, for instance that with w > 

,„uAn „ ,u.out I o {^u,out\^ /on\ 

(fj = Lp^' + (j^ [Lp_^ ) . (20) 
Indeed, the norm of (3^ is given by both overlaps 

= i(<^-'r*,¥':i'")r = i(¥':i'°"**,¥'-'r)i'- (21) 

This square fixes the mean occupation number of 
phonons (found in the in vacuum) for both positive and 
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negative frequency: h^^ = h-ui = \0lj\'^- In fact, since 
each produced pair contains one quantum of frequency 
LO and its partner of frequency —u, there is no need to 
differentiate between these two occupation numbers. We 
emphasize this simple property because it will be lost 
when the coupling to w-modes is no longer neglected. 

A significant difference due to dispersion is that n^^ = 
for LO > Wmax because there is no bounded mode <^"^^ for 
superluminal dispersion (no v'" for subluminal disper- 
sion). 

2. General case with u — v mixing 

Taking into account the coupling between right and 
left-moving modes, the above situation changes as fol- 
lows. 

For LO > Wmax: onc still has n^^ = 0. However, the 
remaining two modes are scattered by the potential, and 
this is described by a single equation 

,„u.in rp u,out I r> ,^v.out /ooA 
foj' =J-i^Vu > (22) 

where |Tljp -t- = 1. One thus has an elastic scatter- 
ing of u and v modes, without spontaneous pair creation. 

For LO < tj^max, unlikc what was found in Eq. (20), three 
equations are now needed to characterize the Bogoliubov 
transformation between in and out modes 

if J = a^^pj + p^i^ (93_^ j +A^ipj , 

As usual, the coefficients are given by the KG overlap of 
the corresponding (normalized) in and out modes, e.g. 

= — ((^^'°"**, (/j"'™). Using the mode orthonormality, 
their normalization immediately follows, e.g. for the first 
equation one gets \a^\'^ — + \Ai^\'^ = 1. In these 

expressions, the minus signs come from complex conju- 
gated, negative norm, modes. 

When working in the in vacuum, the occupation num- 
bers of out quanta, are respectively 

f^LJ I /^CJ I , 

= + = + (24) 

To obtain these expressions one needs to compute the 
in — in expectation value of the corresponding out oc- 
cupation number operator. The three expressions follow 
when decomposing the out operators in terms of in ones, 
and using relations amongst Bogoliubov coefficients. 

Several remarks should be made. First, to get these 
expressions we have treated the quanta emitted to the 
right, which correspond to the outgoing Hawking radi- 
ation and which are described by ip"'°"^, on the same 
footing as those emitted to the left which are described 



by Lp^'°^^ and t/'"'"" . When v is asymptotically constant 
on both sides, there is no reason to treat them differently. 

Secondly, because of the u — v coupling, the numbers 
of u quanta in general differ: n_(^ n^. In fact the 
meaning of n_;j = + is clear. It tells us that 
two channels exist to spontaneously produce u-quanta of 
frequency — w: either through the usual channel where 
the partner is a Hawking quantum reaching x = oo, or 
through the new channel where the partner is a u mode. 
When this B channel is negligible, i.e. when \BJ\^ <C 
one recovers the former situation where 

The B channel was first described in [10], where it was 
also claimed that "this particle creation has absolutely 
nothing to do with black holes." We do not agree be- 
cause, in our settings, there is no reason either to treat 
differently the production of u — w pairs weighted by B^^ 
from Hawking radiation, i.e. the production of u—u pairs 
weighted by (i^^. Indeed, both and n"^ are generically 
non zero when (and only when) the negative frequency u- 
modes Lp^^ exist. ^ More precisely, unless one deals with 
a 2D massless relativistic field, the coefficients B^^ do not 
vanish. In fact, one generically has 7^ even when 
using the dispersionless mode equation of phonons in a 
Eulerian fluid [9, 18, 21]. 

Thirdly, there is a linear scattering between u and v 
modes of positive lo which is characterized by the co- 
efficients and for w < Wmax and by R^^T^ for 
LO > Wmax- This is reminiscent of the "grey body fac- 
tors" which are found when considering four dimensional 
black holes. However the latter are computed on one side 
only of the horizon, whereas here the coefficients relate 
asymptotic modes defined on both sides. 

Finally, we emphasize the fact that Eq. (23) describes 
the scattering whenever a stationary background con- 
tains a single horizon surrounded by two asymptotic ho- 
mogeneous regions. It applies indeed to all dispersive 
mode equations, both for sub- and superluminal rela- 
tions. It is only when the u—v mixing identically vanishes 
that Eq. (23) reduces to Eq. (20). 



E. Wave-packet propagation 

It is useful to represent the above in modes. The de- 
scription of out modes follows without difficulty. Let us 
start with (f^^ in the superluminal case. This mode has 
the following space-time stucture, see Fig. 5. Initially, 



^ It is appropriate to raise the question whether the existence of 
the negative frequency modes implies that the background 
is similar to that of a "black hole" . In this respect it should first 
be pointed out that the definition of a black hole is inherently 
ambiguous in the presence of dispersion. Nevertheless, what is 
always true [in a metric as in Eq. (1)] is that whenever ^p^^ 
exists, its group velocity dQ/dk is smaller than v in some region 
of space, as it is the case in the inside region of a black hole when 
using relativistic fields. 



Figure 5: Schematic space-time behavior of a wavepacket 
made out of ip^^ modes, for superluminal dispersion. 



persion. 



one only has the incoming branch which possesses a unit 
norm, and whose initial wave vector is fc", the largest 
additional real root when the flow is supersonic. At late 
time, one has three branches: the "reflected" negative 
frequency mode described by (/'-w"* whose amplitude is 
a_aj and whose wave vector is fc" , the smallest new root; 
the produced v mode described by (95^'°"*)*, and the pro- 
duced u mode described by (</?^'°"*)*- We leave the de- 
scription of the other in modes and of the out modes up 
to the reader. 

From the above we can already conclude that the emis- 
sion of radiation stops for u) > Wmax- Indeed, in the ab- 
sence of negative frequency partners, both n^^ and nl, 
identically vanish above Wmax- The same is not true 
for the subluminal case. It is thus worth describing 
the in mode v'-i" ™ case, see Fig. 6. Initially, 
one still has only the incoming branch which possesses 
a unit norm and which has the largest positive wave 
vector fc". At late time, one has again three branches: 
the transmitted negative frequency mode described by 
^«,o«t ^j^Qgg amplitude is a-^j, the produced v mode de- 
scribed by and the produced u mode described 
by When w > Wmax, the production of posi- 
tive frequency u modes vanishes {n,^ = 0) but the 'new' 
channel is still open. Thus, in the case of subluminal 
dispersion, one exactly has fi-^ = tV^ for u> > Wmax, 
and no longer = as for superluminal 
dispersion. 

From a conceptual point of view, this completes the 
analysis of the black hole case. What remains to be done 
is to compute the coefficients of Eqs. (23). If one is only 
interested in the flux in the in vacuum, the knowledge of 
the norm of P^i and B^^ is sufficient as they fix the three 
occupation numbers of Eq. (24). 

Before performing their calculation, it is interesting to 
consider flows that engender white holes (WH). Given 



that —v{x) describes a WH geometry when v{x) de- 
scribed a BH geometry, the fluxes emitted by a WH can 
be algebraically related to those of the corresponding BH. 



F. White holes 

We first notice that the mode equation (5) is left un- 
changed under the double replacement v ^ and 
t — > —t. Then, in Fourier transform, the left-moving 
positive norm modes with fixed lu in the WH geometry 
are the complex conjugate of the right movers found in 
the corresponding BH geometry. This means that each 
mode has a group velocity Vg = duj/dk opposite that of 
the corresponding BH mode, so that the BH in modes 
become out modes and vice versa. 

When working in the in vacuum (defined in the same 
way as in the former subsections), and when neglecting 
the coupling to the w-modes, this implies no change in 
the occupation numbers: 



-WH 



(25) 



where is given by |/3ajp obeying Eq. (21). 

Before discussing the physical consequences of this re- 
sult, let us see how it gets modified when taking into 
account the u — v mixing. In this case the occupation 
numbers differ. To obtain them, it suffices to write the 
inverse Bogoliubov transformation of Eq. (23) , and relate 
the coefficients of the inverse transformation to those of 
Eq. (23). One finds that the occupation numbers in the 
WH geometry governed by —v{x) are exactly given by 



"'(aJ 

= It WH 



-i^ I I 



WH 



1/3- 

'^ — OJ'i 



(26) 
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where the Bogohubov coefficients arc computed in the 
BH geometry described by v{x) and where the super- 
script u indicates that l-B^^I^ corresponds to the occupa- 
tion number of the right movers in the WH geometry. 

The equality n^i^ = n-^ expresses that the (total) 
pair creation rate is unchanged, as it must be the case 
because it governs the norm of the overlap between the 
in and out vacua. The inequality ^ means that 
their internal repartition docs not coincide. However, 
when \B^\^ < one has n^H ^ ^wh ^ 

this case, the WH fluxes are essentially the same as those 
of the corresponding BH fluxes. 

Besides this close agreement, the important conse- 
quence of Eq. (26) is that dispersive theories predict that 
WH emit fluxes with well-defined asymptotic properties 
(when starting with the asymptotic in vacuum, which 
is well-defined when v{x) is asymptotically constant). 
When using relativistic theories instead, one obtains an 
endless focusing of null geodesies toward the (past) hori- 
zon which prevents one from getting any outgoing ra- 
diation. (It should be nevertheless noticed that in the 
equilibrium state, in the so called Hartle-Hawking vac- 
uum, the stress tensor of a relativistic field is regular and 
the spectrum of the particles "emitted" by the past hori- 
zon is thermal with the usual temperature, in virtue of 
the stationarity.) 

In spite of the near equality of occupation numbers 
in Eq. (26), there is a major physical difference between 
the WH and the BH cases, which furthermore explains 
why dispersionless and dispersive theories behave so dif- 
ferently. In the WH geometry, the final values of the 
comoving ("proper") frequencies Cl are of the order of 
the UV scale A, whereas, in the BH case, they are of 
the order of uj ^ Th- This directly follows from the fact 
that the values of and k of WH quanta are those of 
the ancestors of Hawking quanta. Because of this when 
taking the limit A — !■ oo, i.e. the dispersionless limit, 
the BH fluxes are asymptotically unchanged whereas the 
WH ones become singular (and ill-defined). This shows 
once more [32] that the relativistic (dispersionless) case 
constitutes an isolated and unstable far as UV 

properties are concerned. 

In brief, we have shown that the ill-definedness of WH 
fluxes in relativistic theories, which is due to the un- 
bounded character of the blue shifting effect, is properly 
regularized by the use of dispersion. (In this we disagree 
with the conclusion reached in [33].) 



G. CJ modes 

We now address the question of the numerical analy- 
sis of the solutions of Eq. (7). To this end, one needs 
to introduce yet another type of mode. We shall call 
them the "CJ" modes, since their usefulness for numeri- 
cal analysis was first recognized in [10]. The basic reason 
to introduce these modes has to do with the control of the 
growing modes. Whereas conceptually we have demon- 




Figure 7: Space-time behavior of v?J = tpw'""': for subluminal 
dispersion. 

strated that the growing modes should be discarded from 
the field operator, when performing a numerical integra- 
tion of Eq. (7), they will be systematically generated and 
will thus completely mask the physical modes. To avoid 
this nuisance, one should use the CJ modes. Even though 
they are defined essentially in the same way for sub- and 
superluminal dispersion, we analyze them separately be- 
cause the information they carry differs. 

1. Subluminal dispersion 

Consider the integration from left to right of Eq. (7) 
starting deep in the supersonic region, i.e. on the left 
of the horizon, where the complcx-fc modes live. When 
reaching the horizon region, the solution is completely 
dominated by the mode that grows when going toward 
the horizon (in our terminology, we called it the decay- 
ing mode, see the Appendix). Therefore, one is effectively 
left with a well-defined mode, up to an overall normaliza- 
tion. The CJ mode is the mode which has a unit positive 
KG norm, and such that only the coefficient of the decay- 
ing mode is non zero in the asymptotic region x — oo. 
For X oo it contains four oscillatory asymptotic modes, 
see Fig. 7. It should be noticed that this is only true for 
quartic dispersion. Indeed for higher order dispersion, 
see Eq. (8), there will still he p—1 growing modes which 
complicate obtaining a stable numerical analysis. 

By direct inspection, one verifies that the CJ mode cor- 
responds to the out mode Using the conventions 
of Eqs. (23), one has indeed 

'P>^='PJ ="a;^a; - Pu, -u; ) +K'PJ -(27) 

Thus one can read off the occupation number of Hawking 
quanta = |/3aiP- However, using the CJ mode, one has 
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We now have to establish the hnk between these coef- 
ficients and the occupation numbers of out quanta in the 
in vacuum. Using the above equations, one can deduce 



and 



\Bj' = \b: 



C\2 



(30) 



(31) 



In the sequel we shall make the approximation (already 
used in [5]) n^^ = |/3^P to compute the modifications in- 
duced by superluminal dispersion. This approximation 
shall be validated by showing that l^ujp <C 1 is veri- 
fied in the whole region of the parameter space where we 
could compute it, and that \Au\^ is decreasing toward 
the inaccessible region. 



Figure 8: Space-time behavior of tp-^, for superluminal dis- 
persion. 

no access to the number oiv quanta = l-B^^p. To have 
access to we shall use other modes, when the growth 
toward the horizon is not too strong. ^ 



2. Superluminal dispersion 

In the superluminal case, the situation is symmetrical 
with respect to the horizon, but also more complicated. 
The diSiculty arises from the fact that the C J mode does 
not correspond to an in nor an out mode, since it pos- 
sesses two incoming and two outgoing waves as shown in 
Fig. 8. Nevertheless it is still defined as the unit positive 
norm solution that is purely decaying (growing toward 
the horizon) in the region where complex-A; modes live. 
To have a positive norm, it must have a negative fre- 
quency. Hence we shall denote it <p^^^■ It possesses four 
branches for x — > —oo with well-defined asymptotic prop- 
erties: a wave packet made out of CJ modes has, at early 
times, the following behavior 

At late times, it behaves instead as 

V^^=l^^^-t+BS{^l'""'y . (29) 
The coefiicients obey |a?^|2-|/3^|2 = \'y^J'^-\B^\'^ = 1. 



^ For high values of oJmax/K the growth of the CJ mode is so strong 
that we were not able to extract information about the other 
modes. Instead, for tJmax/K J; 2, tlic growth is sufficiently mild 
that we could compute accurately the two oscillatory solutions 
and thus extract the complete Bogoliubov transformation. 



IV. NUMERICAL PROCEDURE 
A. Initial conditions 

As explained above, the presence of complex roots im- 
poses the direction of integration and the way the initial 
conditions are fixed. Indeed, when integrating Eq. (7) 
numerically, the mode growing in the direction of integra- 
tion will dominate at some point the oscillatory modes. 
In the subluminal case, to get the CJ mode one must 
integrate from left to right starting from deep inside the 
supersonic region. For the same reason, to get the CJ 
mode in the superluminal case, the integration must be 
performed from right to left. 

In addition, the imaginary part of characterizing 
the CJ mode is gcnerically large (in units of k) even for 
moderate values of D and A = A/k. For instance, for D = 
0.5, A = 50, and lo = k, and for subluminal dispersion, 
|ImA:'^|/K ~ 55. This constrains how far from the horizon 
the initial conditions can actually be set. 

In practice, the initial conditions are thus fixed as fol- 
lows. At a point Xi deep into the supersonic (subsonic for 
superluminal dispersion) region, so that v{xi) is equal to 
its asymptotic value up to the machine precision (10~^^ 
with the C double precision type), the value of the mode, 
its first, second and third derivatives are imposed using 
the fact that in this asymptotic region it is equal to 

=woe^'=^("""'^ (32) 

Mo is taken very small so that the initial exponential 
growth of the mode does not cause an overfiow (the C 
double precision type is limited to the range 10~^°^ — 
10^°^). Since uq cannot be smaller than 10"^°^, \xi\ 
is limited by the imaginary part of fc^. This limita- 
tion causes no problem since, as we saw, the bigger the 
growth of the CJ mode, the less important the preci- 
sion of the initial conditions, the growth itself ensuring 
that any trace of unwanted excitation of the other modes 
quickly becomes smaller than the numerical noise. 
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In brief we have two regimes, either |Imfc'-^|/K 1 in 
which case the growth is strong and we have only access 
to the CJ mode, or \lmk'-'\/n < 1 and the growth is 
weak and we can compute all modes, see Footnote 3. 
Luckily, in the first case, the leading deviations wrt to 
the standard flux can be deduced from the CJ mode. 



B. Extraction of the asymptotic coefficients 

The integration is carried out from Xi to some Xf on 
the other side of the horizon, by an embedded 8th order 
Runge-Kutta-Prince-Dormand algorithm, with a relative 
precision of 10~^^. The real and imaginary part of the 
numerical solution to Eq. (7) thus obtained are stored 
between some Xw and chosen so that it sits in 

the region where v{x) has reached its asymptotic value. It 
is typically taken equal to — Xj. Xf is then chosen so that 
\xf — Xw\ is equal to the period of the component with 
the smallest \k\. which on Fig. 2 is seen to be , the wave 
vector of the f-mode. There is thus no limitation on the 
accessible domain of frequencies to, besides the fact that 
for very low uj, one needs a lot of memory to store (enough 
data to keep track of the component with the smallest 
wavelength, of order A~^. For strongly dispersive cases, 
we could reliably access frequencies as low as 10~^Tfj-, 
where Th — k/2it is the Hawking temperature. 

On the interval [xw,Xf], the mode is a sum of the four 
asymptotic solutions: 



Bogoliubov coefficients are thus given by 



(33) 



Knowing the four roots fcj(w), we extract the coefficients 
Cj by a least-square fitting procedure, tolerating relative 
errors on the various coefficients of less than 10~^. In 
practice the relative precision was usually much better, 
typically 10~^. We are nevertheless limited by the preci- 
sion of the numerical integration: the fit fails when the 
smallest coefficient is smaller than about 10"'^'' times the 
largest coefficient, which means that the smallest compo- 
nent in Eq. (33) is at the level of the numerical error. This 
forbids one to approach arbitrarily close to Wmax since, as 
we shall see, the occupation number quickly drops when 
Lu tiJinax- lu the robust regime, this sets an absolute 
higher bound for the explorable values of w/k. Indeed 
in this regime, f3u/aoj =i e"^'^'^/'*, so P^j = W~^'^auj is 
reached for lj/k ~ 5, i.e. uj/Th — 30. 

The numerical solution to the wave equation is not nor- 
malized, and neither are the asymptotic solutions e''^^'^. 
It does not matter since one only needs the relative norm 
of each component. Taking into account the normaliza- 
tion of modes when v is constant, see Eq. (14), we define 



V 



(34) 



for the positive norm, low momentum u component char- 
acterized by the root k^{uj). In the subluminal case, the 
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|c!i|2 

^2 




\A |2 - 


^2 





(35) 
(36) 



The same equations hold for superluminal dispersion, 
with replaced by (3^ , and A^^ by , and with the 
normalization ijsup given by 

2 I u,out\2 I nt/iu,out\ \ 

Vsup = Ic-'c I ^{kj^ ) 37:^1:^ ' - 



dhj 



(37) 



since the CJ mode is given by Eq. (29) at late times. 

V. RESULTS FOR SUBLUMINAL DISPERSION 
A. General properties of the spectra 

1. Asymptotic energy flux 

The energy flux per du) and in units of k, emitted to 
the right far from the horizon is equal to (denoting by F 
the total energy flux) 



(38) 



The factor 47r^ is added for convenience, so that 1 
for a; — * in the standard case since Th = K/27r. This 
flux is represented in Fig. 9 for A = A/k fixed to 50, for 
values of D from 0.02 to 0.2 (left plot), and for D fixed 
to 0.1 and A varying from 5 to 130 (right plot). 

All curves exhibit a similar shape. When w ^ 0, /a, 
reaches an asymptotic constant value, that differs from 1, 
the dispersionless value, but that becomes very close to it 
whenever Wmax ^ When Wmax/^ is comparable to or 
smaller than 1, the asymptotic energy flux can be either 
smaller or greater than the standard one. In addition, its 
variation with D at fixed A, or with A at fixed D, is not 
monotonic. When w approaches Wmax, the energy flux 
quickly drops to zero, in agreement with the theoretical 
prediction that it vanishes for w > Wmax- 



2. Mixing between right- and left-moving modes 

In Fig. 10, the norm | p that measures the amount of 

elastic scattering between right-moving and left-moving 
modes, is represented as a function of w. The parameters 
for the left and right plots are the same as in Fig. 9. The 
curves display a characteristic shape, with a plateau at 
frequencies much smaller than Wmax, a smooth increase 
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Figure 9: Logarithm of the energy flux density received far from the horizon. The numbers at the top of the plots give the 
value of Wmax/fv for each curve. Left plot: A fixed to 50. D = 0.02, 0.1, 0.2, from left to right. Right plot: D fixed to 0.1. Prom 
left to right, A = 5, 50, 130. In both plots, the squares lie along the thermal energy flux with temperature Th = k/2'k. 
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Figure 10: Logarithm of as a function of uj. In the left plot, A is fixed to 50, and in the right plot D is fixed to 0.1. The 

values of the parameters are identical to those in Fig. 9 



as u) nears Wmax, that gets more pronounced as Wmax 
increases. Then, as for the energy flux, there is a sudden 
fall-off when reaching Wmax- 

These properties can be quahtatively understood by 
considering tlie propagation backwards in time of the 
wave-packet of Fig. 7. As long as the evohition is adi- 
abatic, there is no scattering into left-movers. The scat- 
tering must occur in the relatively well localized region 
of space near the turning point, where the WKB approx- 
imation breaks down. Its position depends on uj: when 



u) <C Wmax, it is very close to the horizon, and as u) in- 
creases, it moves away from the horizon. When it is lo- 
cated in the region where v{x) exits the linear regime 
with slope K, reaches its maximal value. When 
Lo ujmax, the turning point enters into the flat region, 
and goes to zero, as expected. 

In the right plot of Fig. 10, the height of the 
low-frequency plateau significantly grows when A gets 
smaller. It is therefore interesting to further explore the 
behavior of for lower values of A. In Fig. 11, this 
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Figure 11: \A^\'^ as a function of Wmax, for a; = 0.1 x Wmax, for 
D equal to 0.05 (solid line) and 0.1 (dashed line). The range 
in cJmax/K shown corresponds to values of A between 0.5 and 
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Figure 12: Upper plot: logarithm of the energy flux density. 
Lower plot: logarithm of jAi^l'^. D = 0.1 and A — 140 in both 
plots. 



height is shown as a function of the cut-ofF frequency 
Wmax = A/_(_D). The growth seen in Fig. 10 saturates 
when Wmax ~ 0.1 K for both values of D. We also notice 
that l^wP remains smaller than 10~^. 



3. Nonadiahatic effects 

The power n appearing in Eq. (2) cxmtrols the sharp- 
ness of the transition between the linearly increasing flow 
velocity near the horizon and the constant velocity far 
from it. It thus governs the adiabaticity of the wave 
propagation at the transition. The qualitative reasoning 
above suggests that an increase of n should lead to an 
enhancement of and short before Wmax- Fig. 12 

shows that it is indeed the case. 



From the upper plot in Fig. 12, we see that the modi- 
fied /oj contains a small contribution (< 3% for n < 20) 
that shows up when the thermal, exponentially decreas- 
ing part of /a, is sufficiently suppressed [10]. Moreover, we 
see that this contribution is directly related to the nona- 
diabaticity since it increases with n. Instead, u! ^ K/27r 
the low-frequency part of the flux is only slightly affected, 
as is the rapid fall-off for ui Wmax- For I^^^P, a sim- 
ilar enhancement is observed, but it affects all values of 
CO. This is not surprising, since A^^ is entirely due to 
nonadiahatic effects. 



4- Detailed Analysis 

In the next sections we successively want to 

• identify the region of the parameter space {D, A) 
where the outgoing flux /^^ is robust, i.e., where its 
value at w = Th, differs little from the standard 
value = (e — 1)~^. 

• determine, in the robust regime, the behavior of 

both the leading corrections to the thermal flux, 
and of the coefficient lA^^p. 

• analyze the spectral properties away from the ro- 
bust regime. 

• analyze n" when the growth of the C J mode is mild. 

• investigate the properties of the fall-off near Wmax- 

• analyze the integrated energy flux. 



B. Robust regime 

We define the robustness of the outgoing radiation by 
the fact that the energy flux, Eq. (38), differs little when 
evaluated around uj = Th, from the thermal flux ob- 
tained without dispersion: 



f" = 



CO 



Th e-yc^iuj/Tn) - 1 ' 



(39) 



Our aim is twofold. First we wish to characterize the 

region of the parameter space where of Eq. (38) hardly 
differs from . Then we want to determine the scaling 
properties of the modifications. 

To quantitatively study these aspects, we define A// 
to be the relative difference between the modified and 
standard energy flux at w = Th'- 



A 



H 



Jul Ju) 



H 



(40) 



It should be clear that the eriterium Ah <C 1 is local 
in that it is only concerned with what happens for fre- 
quencies near the Hawking temperature. The smallness 
of Ah does not imply that (/ — ,f^)/ remains small 
for all OJ, as can be seen from the upper plot in Fig. 12. 
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Figure 13: A/f of Eq. (40), as a function of Wmax/t for A = 8 
and A = 160. In spite of this factor 20, the deviations are 
extremely similar for all values of Wmax- 
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Figure 14: log|Aj/| and log|Atj=TffP as a function of 
logajmax/K for several values of D. When ct^max > 6, both 
these quantities become linear, with a slope close to —4. 



1- Wmax//t is the most relevant parameter 

In Fig. 13, Ah is represented as a function of Wmax/t 
for A = 8 and A = 160. Although the two values of A 

differ by a factor of 20, both curves stay very close to 
each other even outside the robust regime. This confirms 
that A/k does not govern the robustness of the radiation. 
A is nevertheless relevant since cjmax = A/(D). Hence 
A/k ^ 1 is not a sufficient condition for the robustness, 

but Wmax/ K » 1 is. 

With more details, starting from the right of Fig. 13 
where Wmax/w ^ 1, we find as expected that Ah asymp- 
totically vanishes. It stays much smaller than 1 up to 
t^max/K of the order of 2. For smaller values of Wmax/w, 
Ah becomes positive and reaches its maximal value, 
^ 25%, for LOmax/ K ~ 0.6. Finally, when aj^iax/ k 1/2tt 
(and thus Wmax u) = Th), Ah — » —1, which means 
that the flux vanishes, as already discussed. 

2. Leading corrections in the robust regime 

To determine the scaling properties of the corrections, 
we have represented in Fig. 14, log | A^l as a function of 
logcjmax/'t for several fixed values of D, That is, for each 
curve, A varies along the curve. The minimum value of 
Wmax is such that Ah is already in the region where it is 
small and negative (see Fig. 13). 

We see that, whatever the value of D, for Wmax > 6k 
(logcJinax/K ^ 0.75), log \Ah\ is a monotonic linearly de- 
creasing function of logWmax- The fact that, for different 
values of D, the linear regime starts almost at the same 
Wmax is a further illustration of its relevance. Had we cho- 
sen log A instead, the linear regime would start at very 
different abscissa for each value of D. 

In the linear regime, the slope is of the order of —4, but 
slightly varies with D. The precise values of the slopes are 



-4.42, -4.18, -4.10 for D = 0.2, 0.5, 0.9, respectively, 
so the slope gets closer to —4 for a higher D. Supris- 
ingly, a similar behavior is observed for log|At^=THp- 
The slope is slightly less sensitive to D, as it is equal 
to -4.05, -4.03, -4.01 respectively. 

The fact that the robust regime, Ah <C 1, is already 
reached for Wmax/^ > 2, together with the fact that the 
slope of log Aff is of the order of —4 (and not —1 or —2 
as one might have expected, and as was found for the 
power spectrum in inflationary cosmology [29]) demon- 
strate in a precise and quantitative manner that the ther- 
mal properties of Hawking radiation are hardly affected 
hy UV dispersion. 



3. Role of the asymptotic velocity v+ = —1 + D. 

Figure 15 shows the same quantities, but now with a 
flxed value of A for each curve, while D varies along the 
curves. Linear regimes starting at Wmax/ k — 6 are again 
obtained but the slope is of the order of —3 for log \Ah\, 
whereas it is only —0.2 for log \A^=Th (More precisely, 
the slope for log|A//| is —2.95 for A = 50 and —3.05 
for A = 100, whereas, for log \ A^^Th l"^ ^ it is —0.14 and 
—0.29 respectively.) D thus affects very differently the 
pair creation rate and the elastic scattering. This can be 
understood as follows. Since D fixes the size of the near 
horizon region where the gradient of v can be approxi- 
mated by a constant, D governs the amount of redshift 
suffered by right movers from their turning point to in- 
finity. One thus expects that the deviations from ther- 
mality are highly sensitive to D, and this is indeed the 
case since the log- log power is near to —3. On the con- 
trary, the scattering between right and left movers is not 
related to the size of the near horizon region. Therefore 
one does not expect that log|^ij=TffP be strongly de- 
pendent on logD. The numerical result ~ — 0.2 confirms 
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Figure 15: loglAnl and log|A„^Tffp as a function of 

loga;max/K for several values of A. When a;max > 6, both 
these quantities become linear, log |Ajf | has a slope close to 
3 for both values of A, while the slope of log \A^=Th P is very 



Figure 17: Effective temperature To, as a function of w/oimax, 
for D = 0.1 and various values of A. The dots show the fits 
obtained using the ansatz Eq. (43) . 
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Figure 16: log|Air| and log|Aa;=Tj/P as a function of 
logWmax/K for £) = 0.5 and different values of n. 



this expectation. 



4- Conclusion and comparison with former work 



This result differs from what was reported in [10]. In 
that paper, D was fixed to 0.5, and the scahng of A// with 
respect to A was reported to be characterized by a power 
close to 1. However there is no real contradiction as can 
be understood from Fig. 16, where the same quantities 
are represented as in the previous figures, for D = 0.5, 
and several values of n. We see that, as n increases, 
the superimposed oscillations become more pronounced 
and last until higher values of ciJniax/K. In fact we see 
two regimes. When n is smaller than about 4, the linear 
regime with a slope close to —4 remains. Instead, for n = 
5 the nonadiabatic effects dominate and the linear regime 
is lost. In this sense, the low values of n correspond to 
a superrobust regime, with well-defined first corrections 
to thermality and with a clear scaling in the parameters. 
In [10], the scaling was estimated for a kinked velocity 
profile, which corresponds to the limit n — > oo. In this 
limit, the nonadiabatic effects largely dominate and this 
probably explains why the authors missed the scaling as 
(A)-l ■ 

It is also interesting to note that the behavior of 
log|A^;=Tjj|^ is more robust to changes in n since the 
linear regime exists longer, and the slope and intercept 
are not significantly modified. 



Since the slopes obtained in Fig. 14 and Fig. 15 differ. 

Ah does not scale as a power of ti'max/'^ alone. Remem- 
bering that Wmax = f{D)A, we conclude that deep in 
the robust regime, the first corrections to the thermal 
spectrum can be written as: 

A. = -,(AA/.)x(-A_)\(£), (41) 

where ,g(_D, A/k) is positive and varies only slowly with 
the parameters. 



C. Black hole radiation outside the robust regime 

In the strongly dispersive regime, if one can expect 
large deviations, i.e. 0(1), wrt to the standard flux of 
Eq. (39), one has a priori no idea of what the proper- 
ties of these deviations could be, nor how generic they 
are. Moreover, we do not think that they can be com- 
puted analytically. However with our code we were able 
to compute them when o^max/'* ^ 1, and, to our sur- 
prise, we found that thermality is preserved in the low 
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frequency part of the spectrum. In addition, because of 
the smallness of Wmax/'^; the growth of the CJ mode is 
moderate, and this allowed us to compute the observ- 
able which we could not have access to in the robust 
regime. 



1. Right-moving flux riu; 

To characterize the modified properties of the flux we 
found a posteriori that it is convenient to use the effective 
temperature defined by: 



= = 



1 



(42) 



In Fig. 17, Ti^/Th is represented as a function of 
w/wmax for D = 0.1 and different values of A. Besides 
the expected quick drop of the power when lu — > Wmax, 
two remarkable features emerge: first, for all values of 
"-■'max, is nearly constant for lo < cjniax/lO. Secondly, 
the asymptotic temperature Tq = T^^^o strongly differs 
from Th when uJnia^/K is small. 

These two properties siiggest studying two types of cor- 
rections: a (global) temperature shift, characterized by 
To — Th, and a deviation from thermality, characterized 
by the running of in u). To sort out these two effects, 
we fit the curves with 



Tu;/Th = To/Th + R 



(43) 



The fits are shown as dotted lines in Fig. 17. They coin- 
cide perfectly with the curves at low frequencies and start 
to depart only around /2. The fit parameters 

are found to 
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As expected, when LUrrmx/ is large, Tq becomes very close 
to Th- In all cases, the power s is close to 2, and the 
coefficient R is much smaller than 1. 

We are now in a position to determine to what extent 
the modified spectrum can be considered thermal. To 
this end, let us compare the temperature shift 



Ao 



Tn 



Th 



with the running of T^^, that we characterize by 

To - T^ 



Ato 



(44) 



(45) 



{T^=To is the value given by the true curve, not the fit.) 
We find: 




0)/ K 



Figure 18: Effective temperature To, as a function of w/k, for 
Wmax/w = 0.8. Upper plot: n fixed to 2, A varies; lower plot: 
A = 10, n varies. 
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When Wmax/K is large, these numbers make more precise 
the statement made before that wc recover a Planck spec- 
trum with the standard temperature, truncated when 

— * Wmax- Indeed, both the temperature shift and the 
running in oj are extremely small. 

When Wmax/^ becomes small enough, as it is the case 
for values A < 50 in the table, the departure from the 
standard fiux becomes significant. Nevertheless, there 
is a clear ordering between the temperature shift and 
the running, the former being always much bigger than 
the latter. Thus, as a first approximation, the modified 
spectrum can be seen as a Planckian spectrum with a 
nonstandard temperature, and truncated for lu Wmax; 
and this even far away from the robust regime. 

To complete our analysis, in Fig. 18, T^^ is shown for 
several couples (A, D), such that ciJniax/t remain fixed to 
0.8. To slightly varies (^ 3%) from one curve to another, 
and is thus not controlled only by Wmax- However, the rel- 
ative variations of To are much smaller than the changes 
in A. This shows what we had announced, namely that 
the modification of the spectra is essentially governed 
by i^max/i^ and almost degenerate along the contours of 
constant Wmax/^ in the (A, I?) plane. In fact, when A 
becomes high, Tq saturates at a maximum value (^ Th) 
which is truly fixed by Wmax/ k alone. 

In the lower plot of Fig. 18, the effect of n on T^^ is 
investigated, for Wniax/^ = 0.8 and A = 10. The higher 
n, the higher Tq. The value of To saturates at a maximum 
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Figure 19: Loft plot: as a function of cj/t^max for various 
values of Ljmax/'t, and with D = 0.1 and n — 2. Right plot: 
Tl^^ as a function of uj/k for various values of n, and with 
Wmax/Zt = 0.85, D = 0.1. 



Figure 20: Left plot: relative difference between and 
as a function of a;/ci;max, for cumax/K = 0.085 and 0.85, and 
D = 0.1. Right plot: nZ and function of uj/k for 

Wmax/K = 0.85 and D = 0.1. 



value when n becomes high. We have also checked that 
the dependence of To on n is smaller when A is higher, 
as could be expected by looking at Fig. 12. 

To conclude this subsection, wc point out that the fact 
that the temperature shift of Eq. (44) is much larger 
than the running governed by At^ of Eq. (45) for the 
entire class of metrics we considered severely limits the 
possibility, raised in [27], that "The a;-dcpendcnce of 
the Hawking temperature can be explained by the fact 
that high-energy wave-packets have a different group ve- 
locity than those at low energy and hence the various 
modes 'see' different horizons and thus other values for 
the surface gravity" . Were this picture valid, the shift Aq 
would vanish since low-frequency wave-packets travel at 
the speed c = 1 and thus see the standard surface gravity 
K = dxv\x=o. 



2. Left-moving flux nZ 

When LUmax/K < 1.5, the CJ mode grows only slowly 
toward the horizon. It is thus possible to compute the 
other modes, see Footnote 3. With three independent 
modes at hand, we can compute the nine coefficients of 
the Bogoliubov transformation. Here, we present (for the 
first time in the literature) quantitative results concern- 
ing the production of left-moving particles, governed by 
= \Buj\'^. It is represented in the left panel of Fig. 19 
as a function of w/wmax for D fixed to 0.1 and for three 
values of Wmax/^, namely 0.085, 0.17 and 0.85. We note 
that is, to a good approximation, linear in uj/umax- 
Hence it does not go to zero when w Wmax, in agree- 
ment with the fact that both <^"'°"* and <^^'°"* , and thus 
the coefficient in Eq. (23), remain well-defined above 
'^max in the subluminal case. 



In the right panel of Figure 19 the influence of the 
nonadiabatic parameter n is investigated. As expected, 
the production of left-movers is more important when n 
is higher. However the dependence on n is weak and 
saturates when n becomes high. 



3. WHvsBH fluxes 

Having access to the full Bogoliubov transformation 
of Eq. (23), we can also compute the diff'erence between 
the fluxes emitted by WH and BH. From the algebraic 
treatment, we can assert (using the unitarity of Eq. (23)) 
that the difference of fluxes obeys 

n^"-n. = -{nl''^"-nl). (46) 

However, to actually compute this difference, a nu- 
merical treatment is required. The relative difference 

{n^^ — nui)/nui is represented in Fig. 20 for D = 0.1 
and uimax/n = 0.085 and 0.85. This difference stays be- 
low 10~^. Moreover it decreases when Wmax/w increases, 
i.e. when going toward the robust regime, as one ex- 
pected since it must vanish when the u-v mixing itself 
vanishes. 

For the larger value of Wmax/ k one can see that the rela- 
tive difference changes sign. This can be understood from 
Eq. (46) and Fig. 19, where we have seen that is pro- 
portional to uj while fil]^^ = |i?a)pj shown in the right 
plot in Fig. 20, is constant at low frequencies. Thus one 
always has > at low enough frequencies. Notice 

also that contrary to the u occupation numbers, 
and fiaj, neither of the n^'s vanishes at Wmax, but they 
become equal, ensuring that at all frequencies Eq. (46) is 
satisfied. 
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Figure 21: log/„ (upper plot) and log|A„p (lower plot) as 
functions of log Au /o^max for A = 50 and several values of 
Wmax/K- The legend applies to both plots. The dotted line 
has a slope exactly equal to 0.5 and fits perfectly the Wmax = 2 
curves. All curves have the same slope in their linear region. 
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Figure 22: Integrated flux in units of k^, as a function of 
Wmax//«- Upper plot: n fixed to 2, several values of D. Lower 
plot: D fixed to 0.2, several values of n. The horizontal line is 

the standard value (= l/487r) of the integrated flux obtained 
using the dispersionless theory. 



D. Near Wmax region 

To complete our study of the spectrum, we investigate 
the quick fall-off near Wmax- We find that, both for the 
energy flux and for the fall-off presents a universal 

behavior (for the class of velocity profiles considered), 
independent of A, D and n. Indeed we find that 

f^ = hp{D,X)x .f^, (47) 

V <^max 

where Auj = cjmax — ^, and similarly for I^^^P, with hp 
replaced by a different function of D and A, hA- 

Fig. 21 shows log /^j and log versus log Aw/wmax, 
near lo = ijJmax- A is fixed to 50 and each curve corre- 
sponds to a different value of tJmax- Both quantities are 
linear when oj is close enough to Wmax, with a slope pre- 
cisely equal to 0.5, independently of the value of LOmax- 
Wc verified that this behavior is also independent of A. 
Finally wc have also verified that the value of n, which 
strongly affects the behavior of the spectrum near Wmax, 
as can be seen from Fig. 12, does not affect Eq. (47). 

E. Integrated energy flux 

In this section, we present the main features of the 
integrated energy flux obtained by integrating over u) the 

differential flux /^j. It is computed as a Riemann sum 
and depends on D, A, and n. The sum is computed from 



ioj K— 10"'^, where the asymptotic value of ]^ is reached 
to a good precision, until w/k = uj^^y^j The step 

vhuj/k is 0.01, so that the relative variation between two 
successive points is small. Under these conditions, the 
Riemann sum gives a very good estimate of the integral. 

The upper plot in Fig. 22 shows the integrated flux 

as a function of cjmax/K for several values of D, and n 
fixed to 2. Like Tq and A//, the integrated fiux con- 
verges rapidly to its standard value for Wmax/K > 2, even 
though the flux density differs from the standard one 
for high frequencies. In other words, the nonlocal ro- 
bustness criterium, that the integrated flux be close to 
standard, agrees with the local one defined in Sec. VB. 
This is because the exponential tail of the standard 
contributes only marginally to the integral, so that the 
cut-off at Wmax in the modified flux does not cause no- 
ticeable differences in the integrated flux. 

The lower plot of Fig. 22 shows that when n increases, 
the amplitude of the oscillations in the integrated flux 
is more important, and they last until higher values of 
t^max- This reflects the fact that the nonthermal part of 
yields a nonnegligible contribution to the integrated 
flux, and can be understood looking at Fig. 12 where the 
enhancement of the high frequency part of the spectrum 
through nonadiabatic effects was shown. 
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Figure 23: The u-v mixing coefficient \A^\^ governing the 
difference between and |/3wp, see Eq. (30), as a function 
of Lo, for various values of Wmax- Left plot: D fixed to 0.1. 
Right plot: A fixed to 50. 




Figure 24: Left plot: relative difference of the energy flux 
densities with superluminal and subluminal dispersion, as a 
function of u>/n. The background geometry is fixed, with 
D = 0.1 and n = 2. Wmax/ft varies from 0.085 to 2.39. Right 
plot: relative difference between and the standard energy 
flux for the same set of parameters. 



VI. RESULTS FOR SUPERLUMINAL 
DISPERSION 



A. General properties of the spectra 



(0 /K=0.1 



Smallness of u-v mixing. 



As can be seen from Eq. (30), is not exactly 

equal to fii^. In the following wc make the approxima- 
tion n^j ~ whose validity requires ^ 1. 
To justify this, we first show that the inequality is in- 
deed verified everywhere we can compute is 
shown in Fig. 23 for different values of Wmax- At fixed 
D, on the left plot, it strongly decreases when Wmax/^ 
increases. For instance, it remains smaller than 10"'' for 
Wmax/K = 1-39. At fixed A, on the right plot, it becomes 
monotonically decreasing only for high enough values of 
D, but stays always smaller than 10~^ for A = 50. These 
results strongly suggest that it must be even smaller in 
the region of the parameter space where it is inaccessi- 
ble to our code. As we shall see in the next subsection, 
the important result is that it is always much smaller 
than both the modifications wrt the standard flux, and 
the relative difference between the super- and subluminal 
fluxes. The approximation — fiw thus induces only 
a negligible error. 

Note also that, unlike what was found in Figs. 10 
and 12, |Aa;P does not vanish for cu Wmax in the super- 
luminal case. This is because (^^'O"* and <^^'°"* remain 
well-defined above Wmax so that, at cOmax, A.^ smoothly 
connects to of Eq. (22). 
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Figure 25: Left plot: relative difference of the energy flux den- 
sities with superluminal and subluminal dispersion. a;max//t 
fixed to 0.1; D varies from 0.02 to 0.3. Right plot: relative 
difference between and the standard energy flux for 
the same set of parameters. 



2. Comparison between super- and subluminal dispersion 



In analogy with Eq. (38), we define the flux 



C\2 



(48) 



To our surprise, we found that is extremely close to 
of Eq. (38) when working with the same background 
geometry [i.e. the same value of D and n) and the same 
value of Wmax- This unexpected similarity shows once 
more that Wmax is the relevant parameter for character- 
izing the departure from the dispersionless case. 
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Although Z^"'' of Eq. (38) and belong to two differ- 
ent models, it does make sense to study their difference 
when one works with the same velocity profile v{x), so as 
to have the same behavior in the dispersionless regime, 
and the same value of Wmax- (Under these conditions, 
the two values of A cannot be exactly the same, but are 
close to one another.) In Fig. 24, the relative difference 
ifS ^ fu^^) / fT'' is represented versus w/wmax for a fixed 
D = 0.1, and a scries of increasing values of Wmax- The 
chosen values of uJma,^/ k, correspond to A = 5,50,140 
respectively for subluminal dispersion, and slightly dif- 
ferent values for superluminal dispersion, tuned to get 
the same Wmax/ The relative differences are of about 
5% when Wmax/^ = 0.085, but are extremely small for 
Wmax/K = 0.85 and Wmax/w = 2.39, except in a thin re- 
gion near Wmax- 

The right plot in Fig. 24 also represents the relative dif- 
ference between and the standard energy flux. When 
the spectrum is not robust, as is the case for Wmax/'v < Ij 
the modifications with respect to the standard spectrum 
are much larger than the differences between the super- 
and subluminal spectra. In this sense, the modifications 
are controlled essentially by the value of LOmax/ and 
seem to depend only marginally on the precise nature of 
the dispersion relation. 

This point can be; made more precise with the help 
of Fig. 25, where the same quantities as in Fig. 24 are 
represented, but this time with Wmax/^ held fixed, equal 
to 0.1, and a varying D. The values of A corresponding 
to D = 0.02,0.05,0.3 are 64.84, 16.36, 1.09 respectively, 
for subluminal dispersion (and of course slightly different 
for superluminal dispersion). From the right plot in that 
figure, it is clear that for this fixed value of Wmax/^, large 
variations of {D, A) cause only small variations in the cor- 
rections to the standard flux. In fact, as expected, the 
greater the scale separation between A and k, the more 
precise the statement that essentially only Wmax/^ mat- 
ters and not the precise nature of the dispersion relation. 
This is illustrated by the left plot, where we see that as 
D decreases (and thus as A increases to maintain Wmax/ n 
fixed), the relative difference between and be- 
comes smaller and smaller. 

The smallness of the differences in the energy flux be- 
tween sub- and superluminal dispersion teaches us several 
lessons. First, as soon as uj^ox/k, > 1, the quantitative 
agreement, to a fraction of a percent, of the superluminal 
results with the subluminal ones (at least for frequencies 
not too close to Wmax) means that the global properties 
for a given Wmax of the spectra for sub- and superluminal 
dispersion, such as the modified temperature, the run- 
ning of the temperature, and the total energy flux, arc 
the same to a very good approximation. Thus, we do not 
repeat in the case of superluminal dispersion the study 
performed in Sec. VC and Sec. VE, since the results arc 
qualitatively and quantitatively the same. Secondly, the 
robustness criterium and the range of robustness are the 
same for both types of dispersion. That is, the radiation 
is robust for Wmax/ k > 2 in both cases. 



B. First corrections to HR in the robust regime 
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Figure 26: Scaling of with Wmax/K, for fixed D (upper 
plot) and for fixed A (lower plot). 

The above agreement between and suggests 
that the scaling of the leading corrections with respect 
to the standard spectrum are also similar for sub- and 
superluminal dispersion. This point is investigated in 
Fig. 26, where the logarithm of 

^% = {f2-f")lfSU=T„ (49) 

is represented versus the logarithm of Wmax/ « for several 
fixed values of D and A. 

Like in the subluminal case, a "superrobust" regime is 
reached for Wmax/^ ^ 6. The slopes differ also slightly 
from the subluminal case. In the upper plot, we find 
-3.46,-3.79,-3.88 for D = 0.2, D ^ 0.5, D = 0.9 re- 
spectively. In the lower plot, the slopes are —2.42 and 
—2.48 for A = 50 and 100 respectively. The scaling in A, 
at fixed _D, of the first corrections is still with a power 
close to —4 in the superluminal case, while the scaling in 
f+(D) at fixed A is with a power between 2 and 3. It is 
also worth mentioning that is negative in the super- 
robust regime, as was A|J^*. This is opposed to the first 
corrections to the inflationary spectra, that have different 
signs for sub- and superluminal dispersion (see [29]). 



C. Production of left-movers 

Figure 27 shows the occupation number for super- 
luminal dispersion for D = 0.1 and the same values of 
Wmax/w as in Fig. 19. As in the subluminal case, it de- 
creases when cijjnax/t increases. It is proportional to lo for 
u> < 0.1 X Wmax, but Contrary to the subluminal case, it 
goes to zero when w approaches Wmax- This is as expected 
since the negative frequency partners of these u-quanta 
exist only for uj < Wmax- 
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Figure 27: as a function of w/wmax, for Wmax/w from 0.085 
to 0.85, with D fixed to 0.1. 



D. Near Wmax region. 

Finally, we have investigated the behavior of the fluxes 
near the maximum frequency. We find that both riP^ and 
/5 decrease like (t^max - ^Y^"^ near Wmax, as and 
\A'f'f did (see Eq. (47)). 



E. Comparison with Ref. [34] 

In Ref. [34], when considering a timc-dcpcndcnt col- 
lapsing background metric, the authors reached the con- 
clusion that Hawking radiation is no longer stationary 
when dealing with a quartic supcrluminal dispersion re- 
lation. We believe this conclusion is not due to the time- 
dependent character of the geometry, but is rather due to 
the use of a geometric optic approximation, in Eq. (4.14), 
in a regime where it is not valid, as indicated by the nu- 
merical results of [21]. However a detailed calculation of 
the fluxes remains to be done. 



VII. CONCLUSIONS 

We have studied the properties of the fluxes emitted by 
acoustic black holes when taking into account the effects 
induced by sub- and supcrluminal dispersion relations. 
We focused on one dimensional stationary flows which 
possess two asymptotic regions. In this case, unlike what 
is found for gravitational black holes, the fluxes emitted 
on both sides of the horizon are well-defined observables. 

At the theoretical level, we showed that the dimension- 
ality of the space of asymptotically bounded stationary 
modes is three, and that this space is complete. This 
guarantees that the 3x3 Bogoliubov transformation of 
Eq. (23) relating in and out modes is unitary. It reduces 
to a 2 X 2 matrix only if left and right-moving modes 
completely decouple, which is a degenerate case, proba- 



bly never found in condensed matter systems. We then 
showed that because of dispersion, there exists a critical 
freqTK!iicy u;inax above which no radiation is emitted, and 
the 3x3 transformation simplifies into Eq. (22). This 
frequency depends both on the UV scale A characterizing 
the dispersion, and on the asymptotic properties of the 
background geometry. Because of the dispersion as well, 
the fluxes emitted by white holes are well-defined and 
regular. Moreover we established how they are related to 
those of the corresponding black hole, see Eq. (26) and 
Fig. 20. 

We have computed numerically the fluxes in the case 
of quartic sub- and supcrluminal dispersion. The main 
result of our analysis is the following. The deviations wrt 
the standard dispersionlcss fluxes arc highly degenerate 
along the lines of constant frequency Wmaxj see Figs. 13 
and 18, and, as a corollary, are essentially governed by 

'•-'max/ 

Our second result concerns the characterization of the 
range of the robustness, i.e. the range of the parameters 
within which the deviations are small. We first found 
that for Wmax/K > 6, there is a superrobust regime in 
which the leading deviations rapidly decrease as a power 
law in K/wmax, where the power is equal to 3 or 4 depend- 
ing on which parameter is held fixed, see Figs. 15 and 14. 
We also showed that for Wmax/^ > 2, the deviations are 
already smaller than a percent. 

Our third important result concerns the absence of a 
significant running of the effective temperature, Eq. (42), 
even when leaving the robust regime. The low-frequency 
part of the fluxes, i.e. uj < 0.1 x Wmax, is indeed well char- 
acterized by a common temperature that significantly 
differs, for ujra&^/ k < 1, from the standard temperature 
K/27r but that hardly changes with w, see the table after 
Eq. (42). This implies that the effective temperature 
is probably not related in a universal manner to the gradi- 
ent of the flow velocity v evaluated at some a>-dependent 
horizon, in agreement with the remarks made in [10]. 

Fourthly, we have also presented quantitative results 
concerning the production of u-quanta in the domain 
Wmax/w ^ 1- We found numerically that n^, is propor- 
tional to u) at low frequencies, see Figs. 19 and 27. 

Finally we showed in Sec. VI that the deviations ob- 
tained using sub- and supcrluminal dispersion are un- 
expectedly similar when using theories having the same 
t^max- This means that the deviations wrt to the stan- 
dard fluxes obtained with each theory are much larger 
than the difference between these deviations, typically by 
a factor of more than 10, see Fig. 25. However these con- 
clusions are probably related to the fact that our velocity 
profiles are symmetrical wrt the horizon, and possess two 
asymptotic regions. When considering profiles describing 
gravitational black holes, we expect that the deviations 
due to sub- and supcrluminal dispersion will differ more 
significantly. We also expect our results for subluminal 
dispersion to apply to gravitational black holes because 
the CJ modes essentially live in the outside region, as can 
be seen in Fig. 7. 



22 



Acknowledgments 

We would like to thank the organizers and the par- 
ticipants of the workshop "Towards the observation 
of Hawking radiation in condensed matter systems" 
(http://www.uv.es/workshopEHR/) held at IFIC in Va- 
lencia in February 2009, where the extension of this work 
to Bosc-Einstcin condensates was presented, for many in- 
teresting remarks. We are grateful to Jihad Mourad for 
discussions about the asymptotic bounded character of 
physical modes, and to Roberto Balbinot and Ted Ja- 
cobson for useful comments. 



Appendix A: ROLE OF DECAYING MODES 

In this appendix, our goals are, first, to show that in 
the black hole metric of Eq. (1), when using the disper- 
sion relation Eq. (8) and for uj < Wmax, the space of spa- 
tially bounded modes has dimension three, and secondly, 
to understand the physical roles of the decaying modes. 
In Sec. IIIC, when considering homogeneous flows, we 
saw that both growing and decaying modes should be 
discarded from the field operator. When the flow is only 
asymptotically constant, as it is the case in Eq. (1), the 
symmetry between growing and decaying modes is bro- 
ken because a mode that grows toward the horizon (we 
called this a 'decaying' mode) does not necessarily grow 
in the other asymptotic region. 

In the asymptotic supersonic region x — oo of 
Eq. (1), the wave equation reduces to the same form as 
in an homogeneous flow with velocity Thus, in this 
region, the general solution is a combination of all the 
exponential solutions that were present in the homoge- 
neous case: four oscillatory solutions, p — 1 decaying so- 
lutions, and p—l growing solutions. This set defines the 
m,athem,aUcally complete basis of modes (there- 
after called the "left" basis), where a labels the roots of 
Eq. (8) and i stands for "left", such that each 65" has 
only one exponential component e^^'--^ in the asymptotic 
supersonic region. When restricting ourselves to spatially 
bounded modes, the 6f,'" that correspond to growing com- 
ponents must be of course discarded. So one is left with 
p -|- 3 independent modes bounded for x ^ oo, num- 
bered by a = l..p + 3, . 

In the subsonic asymptotic region, in general, these 
modes are not bounded because they are combinations 
of all exponential components, solutions of Eq. (8) with 



velocity w+. On that side, one has p growing and p de- 
caying modes, and only two oscillatory modes. In order 
to get modes that are everywhere bounded, one must 
consider combinations of the form 

p+3 

^i=E^ae- (Al) 

a=l 

The condition that the coefficients of the growing com- 
ponents in the subsonic region be all zero yields p con- 
straints on the p -|- 3 coefficients . Therefore only three 
linearly independent bounded combinations exist. In 
other words, the dimensionality of the space of spatially 
bounded solutions is three for < a; < Wmax- This ex- 
plains why only one extra- mode was added in Eq. (19). 
It is also important to note that the p constraints in no 
way impose that the coefficients of the decaying modes 
be zero, and in fact the bounded modes, and in particu- 
lar the ones constituting the in and out bases, generically 
contain decaying modes. 

The physical role of these decaying modes is to 'dress' 
the in and out modes in the region where v(x) varies, in a 
way similar to the dressing of atoms by a local and non- 
propagating polarization cloud when these are coupled 
to a radiation field [35]. Thus this dressing by decaying 
modes affects the properties of local observables (such as 
correlation functions) in the region where v varies. 

For completeness, a word on the construction of the 
in and out bases is in order. Just as we constructed the 
"left" basis, we can construct a "right" basis {6^'"}, from 
the asymptotic exponential solutions in the subsonic re- 
gion. Looking at the space-time properties of wavepack- 
ets made out of the asymptotically oscillatory modes in 
both and {h[^°'}, one can reclassify them into a 

group of three in modes and one of three out modes (still 
unbounded at this point). Each of these two groups mixes 
"left" and "right" modes. The dimensionality three of the 
space of bounded modes then ensures that we can find 
bounded combinations with only one in (respectively out) 
mode with a nonzero coefficient. These combinations, 
suitably normalized and possibly complex conjugated to 
get a positive norm mode, yield the in (resp. out) ba- 
sis used in the text. Note also that since there are only 
two oscillatory solutions in the subsonic region, there ex- 
ists one (and only one, up to an overall factor) bounded 
mode that is purely decaying in this region. This mode, 
correctly normalized, is the CJ mode, here generalized to 
arbitrary (polynomial) dispersion. 
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